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Oblateness-Perturbed Orbits by Velocity-Correspondence Variations 

V. E. BLEICK 

U. S. Naval Postgraduate School, Monterey, California 

(Received 1966) 

ABSTRACT 

This paper presents an elementary treatment of the first 
order differential effects of the earth's oblateness on a 
close satellite using the simple notion of a varied orbit. 

The usual result of this approach is that the radial vari- 
ation contains a secular term which is unbounded for infi- 
nite time. The standard method of celestial mechanics for 
removing this difficulty is to analyse the perturbed orbit 
as an ellipse whose shape and space orientation are func- 
tions of time. It is shown here that the secular term may 
be avoided more simply by relating points on the varied and 
unvaried orbits by a type of radial velocity correspondence 
instead of the usual time correspondence, and by making 
the varied orbit osculate at a latus rectum chord end point 
of the unvaried orbit. 
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INTRODUCTION 



Let the origin 0 of a non-rotating rectangular coordinate 

frame Oxyz be at the center of the earth with the Oz axis point 

ing toward the North pole. Let the unit mass earth satellite at 

S in Fig. 1 have the coordinates x,y,z, Adopt the equatorial 

radius of the earth as a unit of length. If it be assumed that 

the earth's mass distribution has axial symmetry about the Oz 

polar axis, the gravitational potential V of the satellite may 

1 

be expanded in a series of zonal harmonics. Retention of the 

two lowest order terms in this series gives 

V=-n[ p" 1 +e p" 3 ( 1 -3 z 2 p" 2 ) ] ( 1 ) 

2 2 2 2 

where (i is the gravitational constant and p =x +y +z .A spher- 
ically symmetric distribution of mass corresponds to e=0 in (l) 

2 

Lecar, Sorenson and Eckels have determined that the earth's 

_3 

oblate-spheroidal shape gives rise to the value e=0. 541x10 
squared units of length. The problem is to find the variations 
of an unperturbe'd Keplerian elliptical orbit, corresponding to 
e=0, caused by the small e perturbation term in (l). 

Let z cosl-y sinl=0 be the equation of the unperturbed orbit 
plane Oxv of. Fig. 1 , passing through the Ox axis with inclina- 
tion angle I to the earth's equatorial plane. Introduce the 
cylindrical coordinates r,^,w shown in Fig. 1 with the Ow 
cylinder axis perpendicular to the Oxv plane. The equations 

x=r c o s 5 , ( 2 ) 

y= r cosl sin£-w sinl, 
z=r sinl sin^+w cosl, 
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FIG. 1, Cylindrical coordinates based on the 
Oxv plane of the unperturbed orbit. 



relate the x,y,z and r,?,w coordinates. The kinetic energy T 
and potential energy V of the unit mass perturbed satellite at 
S are then 



and 



T=£(r 2 +r 2 ? 2 +w 2 ) 

13 2 5 

V=-u[p- +ep - -3e(r sinlsin?+w cosl) p - "j 



2 2 2 

where p =r +w . 



(3) 

(4) 



VARIATION OUT OP ORBIT PLANE 



The Lagrangian function T-V may be constructed from (3) and 
(4) , and the equation of satellite motion in the w coordinate 
found to be 



O *N O A P 

w+|rp - w=-ne-g^[ p - -3(r sinlsine;+w cosl)^p - ]. 



(5) 



Using the fact that w=0 in the unperturbed orbit, one finds 
from (5) that the first order time-correspondence variation 6w 
of the perturbed satellite out of the unperturbed orbit plane 
Oxv of Pig. 1 satisfies the equation 

6w+ur - ' ) 6w=-3 enr - sin2Isin5 . (6) 

Note that the operators 6 and d/dt commute in (6) since 6w is 
defined to be a time-correspondence variation; i.e., 6w, w, r 
and § are all defined for the same value of time. Let §=U) be 
the argument of the perigee P of the unperturbed elliptical 
orbit of Pig. 2 measured from the ascending node N on the Ox 
axis. Let §=\=uj— be the argument of the latus rectum chord 
end point L which precedes perigee. If the satellite position 
argument a is measured from L the equation of the unperturbed 
Keplerian orbit may be written as 

r=p/(1+e sina) (7) 
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FIG. 2. Unperturbed Keplerian orbit. 



where e is the orbit eccentricity and p is the semi latus rec- 
tum. The constant angular momentum of the unit mass satellite 



in this orbit is 

h=r^a=(pn)^=2nab/T (8) 

where a,b are the semiprincipal axes and T is the period of the 

orbit. Equations (7) and (8) may be used to replace time by a 

as the independent variable and to replace 6w by r 6w as the 

dependent variable in (6) so that it reduces to 

(r~ 6w) "+r 6w=-3ep sin2I(l+e sina) sin( X+a) (9) 

where the primes indicate differentiation with respect to a. The 

simplicity of presentation of this paper depends on making the 

varied orbit osculate the unperturbed orbit at the latus rectum 

chord end point L. Equation (9) is therefore integrated with the 

initial conditions 6w(0)=6w' (0)=0 at a=0 to obtain 
2 

■^^^j=3acos^-cos\ (3 sina+8esin^^a)-4esin\sinasin^-|a . (10) 

The term proportional to a in ( 10) is a secular term which is 
unbounded for infinite time. To remove this undesirable term 
consider a plane through the origin 0 in Pig. 3, which is init- 
ially coincident with the Oxv plane of the unperturbed orbit at 
a=0, and which precesses around the Oz polar axis to the west 
at the differential angular rate 6cp always maintaining a con- 
stant inclination angle I with the Oxy equatorial plane. The 
equation of this precessing plane is 

-xsinlsin6cp-ysinlcos6cp+zcosl=0 ( 1 1 ) 

where 6cp=0 when a=0. The distance 6w* of the perturbed satellite 
from the precessing plane may be found by substituting its coor- 
dinates x+6x*y+6y , z+6z from (2) into the left member of (11)* 
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FIG. 3. Satellite variations 6 w out of the un- 
perturbed orbit plane and 6w* out of 
the processing perturbed orbit plane. 



noting that w= 0. Using the approximations sin6cps«6cp and cos6cp?«1 , 
and neglecting other second order differentials, one obtains 

6w*=6w-rsinIcos§6cp. (12) 

On comparing ( 1 2 ) with (10) it is seen that the distance 5w* is 
periodic and contains no secular term if 

6cp=3ectcosl/p . (13) 

¥e conclude that the mean plane (ll) of the perturbed orbit 
precesses around the Oz axis at the angular rate 

6<p=3 eacosl/p , (14) 

and that the satellite has periodic excursions 6 w* from this 
plane given by 

6w*=6w-[3earsin2Icos(X+a)/2p^] . (15) 

The relation between 6w*/r and a is shown graphically in Fig. 4 
for a few values of X and e. The trace 0M of the precessing 
plane ( 1 1 ) in the equatorial Oxy plane has the equation y=-x6cp, 
so that the ascending node N of the perturbed orbit in Fig. 3 
moves westward along the equator through the angle 

6cp( 2tt) =6rrecosl/p^ (16) 

on successive transits of the equator. This nodal regression 
from to is shown in Fig. 5. If the node corresponds to 
a on the unperturbed orbit one must not conclude that the next 
ascending node corresponds to <x+2tt since the variations Ar 
and Act in the unperturbed orbit plane of Fig. 3 have not yet 
been considered. It is shown later in (46) that the point Q of 
Fig. 5, corresponding to a+2n, is at a positive or negative dif- 
ferential angle Ai(f ( 2 tt ) of order e from thus justifying the 
assignment of 6cp(2rr) of (16) to the angle N^O^. 
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FIG. 4. Relation between r -1 6v* and a. 




FIG. 5. Regression 6ep(2rr) of the ascending 

node in the equatorial plane and pre- 
cession Ai(r ( 2 tt) of the perturbed orbit. 



VARIATIONS IN THE ORBIT PLANE 



The gravitational force on the satellite has no moment around 
the earth's polar Oz axis because of the axial symmetry of the 
earth's mass distribution assumed in (l). Hence the moment of 
momentum of the perturbed satellite around the Oz axis is con- 
served. The moment of momentum vector is found most easily by 
first computing its components in the Oxw frame of Pig. 1 , and 
then transforming them to the Oxyz frame. It follows from (8) 
that conservation of the Oz component of moment of momentum of 
the unit mass perturbed satellite requires 



• • * 2 * 

[ (vr-rv) cos§-wr§sin§] sinl+r ^cosI=hcosI 



(IT) 



where it is assumed that the orbit osculates the unperturbed 
orbit (7) at §=X , i.e. a=0. On expressing (17) in terms of first 
order time-correspondence variations from the w=w=0 unperturbed 
orbit, and using (8) to replace time by a as the independent 



variable, one obtains 



2r - ^ 6r+6a ' =tan!cos^g^-( zTZZTe ) =Fi 



where, using ( 1 0 ) , 



-p 2 r, 

• 2 t 
esin I 



5 da v rcos§ y 1 



=2e[sin Ja(2cosa+cos2a) sin2X+sin acos2X]+3 sinasin(2X+a) . (19) 



(18) 



From (3) and (4) it is seen that conservation of the total 
energy T+V of the unit mass perturbed satellite requires 

i(r^+r^a^+w^)-n[p“^+ep'"^-3e(rsinIsin§+wcosl)^p'"^] . (20) 

=E-kiep - ^ (l -3sin^Isin^\) 

where it is assumed that the orbit osculates the unperturbed or- 
bit ( 7 ) of total energy E at ?=X , i.e. a=0. On expressing (20) 



in terms of first order time-correspondence variations from the 
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w=w=0 unperturbed orbit one obtains 
r6r+(ra^+(ir ^)6r+r^a6a=(ie[r "^-p ^+3sin^l(p~^sin^X-r - ^sin^?) ] . (21) 
On using (8) to replace time by a as the independent variable 
(21 ) becomes 

r"*^r 1 6r 1 +(r - ^ +p~^ )6r+6a 1 =F^ (22) 

where 

F2=e[l-r^p~^-3sin^I(sin^§-r^p - ^sin^X) ]/pr. ' (23) 

If the rather complicated Eqs . (18) and (22) are integrated as 
they stand it will be found that 6r contains a secular terra, 
exhibited later in (53), which is unbounded for infinite values 

of a. The secular nature of 6r raises doubt about its ability to 

give a sufficiently accurate description of the radial perturba- 
tion when a is more than a few multiples of 2tt, and makes for 

unwieldy mathematics in the determination of 6a. Sterne has re- 
moved the radial secular term in the present problem by a stan- 
dard method of celestial mechanics which analyzes the perturbed 
orbit as an ellipse whose shape and space orientation are set up 
from the beginning to be functions of time. The secular term dif 
ficulty is avoided here in a simpler fashion by relating points 
on the varied and unvaried orbits by a type of radial velocity 
correspondence instead of the usual time correspondence implied 
by the use of the symbols 6r and 6a. We assume that the differ- 
ential change in radial velocity Ar from a point on the unper- 
turbed orbit to a corresponding point on the varied orbit is 

Ar=6r+rAt=rF(a)/r (24) 

where At is the corresponding time difference. The function 
F(a) , of order e, is chosen later to avoid a singularity. The 



function F will also be chosen so that exact radial velocity 
correspondence, Ar=0, is achieved at perigee and apogee, thus 
facilitating a discussion of the precession of the apsidal line 
of the perturbed orbit. The symbol A, used in this section to 
denote the velocity-correspondence variation operator, does not 
commute with the operator d/dt since the correspondence of po- 
ints defined in (24) is not based on equal values of time. The 
velocity-correspondence variations in r, a and a, consistent 



with (24) , are 



Ar=6r+rAt , 


(25) 


Aa=6a+aAt , 


(26) 


and 

Aa=6a+aAt . 


(27) 


Note that At may be found from (25) after Ar-6r 


has been deter- 


mined in (53). Note also that Aw=6w+wAt=6w and 


Aw*=6w* since 


w=0 on the unperturbed orbit. 

Equations (7), (24), (25) and (26) yield 
Ar=6r+[r ' 6r ' / (r 2 p~^ -r) ]+F 


(28) 


and 

Aa=6a+[6r'/( r ^p"" 1-r ) ]+F/r* . 


(29) 


Subtraction of (18) from (22) yields 

r'r“ 2 6r'+(p" 1 -r“ 1 )6r=F 2 -F 1 . 


(30) 


Substitution of (30) into (28) yields 

Ar=F+[(F 2 -F 1 )/(p" 1 -r“' 1 )]. 


(31) 


Elimination of 6r' from (28) and (29) yields 

Aa-6a=(Ar-6r)/r ' . 


(32) 



Diff erentir • ion of (32), together with the use of (7), (18) 
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and (28) , yields 



^Aa=F 1 +(r 2 /r' )[Fr 2 tana+~^( Ar/r 2 ) ] . (33) 

Use of (7), (19) and (23) gives 
• -P ( F 2 -F 1 ) 2 2 2 2 

;_y - sin iLsin acos2\+( sina+2sin §acos acscct) sin2\ ] (34) 

e (p -r ) 9 9 

+(3sin Isin \-1 ) [ 1 +( r/p) +(r/p) 

7T 

which is singular at a=!0. To avoid the same singularity in Ar 



choose F of (31 ) to be 



2 2 
F=esin Isin2\ ( cosa-1 ) cos^a/psina. 

The differential change in radial velocity Ar is then found 

from (24) , with the aid of (7) , (8) and (35) , to be 



(35) 



# * "1 2 

Ar=ectp"" sin Isin2\(l-cosa) cosa, (36) 

giving exact radial velocity correspondence at osculation, peri- 
gee and apogee. The differential change Ar is found from (31), 

(34) and (35) to be 

Ar=ep” { sin I ( sin acos2\+sinasin2\) (37) 

-(1-3sin 2 Isin 2 X) [l +(r/p)+(r/p) 2 ]} , 

giving 

Ar(0)=3ep - ^ (3sin 2 Isin 2 \rl ) . (38) 

The variation Aa(0) at osculation, where 6r(0)=6a(0)=0, is 
found from (7) , (32) and (38) to be 

Aa(0)=3e(l -3 sin 2 Isin 2 x)/p 2 e . (39) 

Equation (39) shows that the notion of a velocity-correspondence 

V 

variation cannot be used for e=0 as one would expect. An inte- 
gration of (33) with the aid of (7), (19), (35) and (37) yields 
Aa-A<x(0)= (40) 

£ep _2 £ [ ( 2 -3 sin 2 l) (3a+4esin 2 ^a) +sin 2 Isina(9cosacos2X-2e - ^ sin2X) ] + 
2sin 2 Isin 2 £a[2e ( cos2a+2cosa-2e~ 2 ) cos2X-(9+5cosa+4esina+esin2a) sin2X]) , 
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periodic in cr except for the secular term 

3ea(2-3sin 2 l)/2p 2 . (41 ) 

The change in a+Aa from a=^rr to a=3rr/2 may be calculated from 
(40) to be 

TT-Hgep _2 [3TT(2-3sin 2 l)-4(2e+e’"^ )sin 2 Isin2X] (42) 

so that the apsidal line of the perturbed orbit does not pass 
through the earth's center unless the last term of (42) happens 
to be zero. The complex variable describing the projection of 
the perturbed satellite orbit in the Oxv plane of Pig. 3 is 

( r+Ar) exp[ i ( X+a+Aa) ] (43) 

if the Ox axis is taken as the axis of the real part. Let OM 
of Fig. 3 be the trace of the precessing perturbed orbit plane 
(11) in the earth's equatorial plane. Let Aijf ( a ) of Fig. 3 be a 
differential angle measured from the trace OM to a moving ref- 
erence radial line OA in the precessing plane. Then the complex 
variable describing the projection of the perturbed satellite 
orbit in the precessing plane relative to the reference radial 
line OA as the axis of real numbers is 

(r+Ar)exp i0 (44) 

where the angle 0 of Fig, 3 is 

0=X+a+Aa+6cpcosI-Ai|i . (45) 

On substituting Aa from (40) and 6cp from (13) into (45) it is 

found that the perturbed orbit projection (44) is periodic re- 
lative to the moving line OA and contains no secular term if 

A\Ji=3ea(4-5sin 2 l)/2p 2 . (46) 

Note that the rate of precession 

(d/dt) A\|/=3ea(4-5sin 2 l)/2p 2 (47) 
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of the periodic perturbed orbit (44 to 46) relative to the trace 



OM is direct or retrograde according as sinl is less than or 
greater than (4/5)^. The arcsin(4/5) 2 is 63°26'. The preces— 
sion Ai|i (2rr)=3eTT(4-5sin I)/p of the periodic perturbed orbit in 
the precessing plane between the times of two successive ascend- 
ing nodes and ^ has been pictured earlier in Fig. 5. 

VARIATION IN NODAL PERIOD 

The change from the period T of the unperturbed orbit to the 
nodal period T+AT between the two successive nodes N^ and ^ of 
Fig. 5 can be calculated only in terms of time-correspondence 
variations. Let At represent the differential time changes from 
the points z=6z at a=-X and a=-X+2n on the varied orbit to the 
corresponding nodes N^ and ^ at z=0 on the varied orbit, so 
that 

6z+zAt=0. (48) 

The differential change in the nodal period is then 

AT=A t ( -X +2tt ) -At ( -X ) (49) 

=[6z(-X)-6z(-X+2tt) ]/z(-x) . 

Substitution of (2) into (49) gives 

a(-x) AT=[ 6 a+r - 1 cotl 6 w]“^ + 2 n . ( 50 ) 

Equation (10) yields 

J 

r - ^ cotl6w| ~^ + 2 TT =- ^ eTT P - ^ C0S ^^ • ( 51 ) 

The calculation of the time-correspondence variation 6a, 
evaded thus far, is now required. Eliminate F 2 ~F^ from (30) 
and (31 ) to obtain 

^-[ (6r-Ar) seca]=-(Ftana+^^-) seca. 

Integrate (52) with the aid of (37)* using (38) and 6r(0)=0 



( 52 ) 



as initial conditions at osculation, to obtain 
(6r-Ar)seca = 

e e (33in 2 I 2 Sin 2 X-1) [ecosa[(4 - e 2 ) s +(1 - e 2 ) ri ] - e( ^ e 2 ) _a (1 _ e 2 ) 2 



( 53 ) 



p(l-e ) 



+3(l-e 2 ) _ 2[a+^ — 

n=1 



oo 



-ep sin z l[4sin z §acos2X+sinasin2\] 

where a and b are coefficients in the Fourier series 
n x n co 

( 1 -e z ) 8 ( 1 +esina) - — g a Q + \ (a^cos na+b n sin na) (54) 

n=1 

and are defined by the complex number 

a n +ib n =2(ie)" n [l-(l-e 2 )^] n . (55) 

Comparison of (37) with (53) shows that the time-correspondence 
variation 6r contains a secular term proportional to acosa, and 
also shows the simplicity of the velocity-correspondence varia- 
tion Ar versus the complicated 5r. Use of (32), (40) and (53) 
gives 

6a | I^ + 2TT = ^ TTe ^ ) ""^^ (3 sin 2 Isin 2 X-1 ) ( 1 -e 2 ) - 5/ 2 +^sin 2 I-1 ] . (56) 

Use of (8), (50), (51) and (56) then shows that the fractional 
change in the nodal period is 

2-r • 2, „ . 

(57) 



• 2i 



a sin na+b (1-cos na) 



n 



n 



■]} 



at 3e r 3sin 2 Isin 2 X-1 ,5. . 2 T ,0 

> L (l _ e 2 } 5/2 + 2 sin I " 2J ‘ 



^ ab(l-esinX) 2 



CONCLUDING REMARKS 



The notion of varying an orbit by velocity-correspondence 

variations instead of time-correspondence variations has been 

shown to lead to a simple discussion of first order oblateness 

effects. The variations obtained are less complicated -than those 

3 

derived by an analysis in terms of an ellipse whose shape and 



1 2 



space orientation are functions of time, and have a more direct 
physical interpretation. The analysis can be carried through for 
a general point of osculation at the price of an increase in 
mathematical complexity. 
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FIGURE CAPTIONS 

FIG. 1 . Cylindrical coordinates based on the Oxv plane of the 
unperturbed orbit. 

FIG. 2. Unperturbed Keplerian orbit. 

FIG. 3. Satellite variations 6w out of the unperturbed 

orbit plane and 6w* out of the precessing perturbed 
orbit plane. 

i 

FIG. 4. Relation between r 6w* and a. 

FIG. 5. Regression 6cp(2rr) of the ascending node in the equa- 
torial plane and precession Ai(f ( 2 tt ) of the perturbed 



orbit . 
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